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IMPINGEMENT OF WATER DROPLETS ON WEDGES AND DOUBLE-WEDGE AIRFOILS 

AT SUPERSONIC SPEEDS ^ 

By John S. Serafini 


SUMMARY 

An analytical solution is j^cesented for the eyiiations of motion 
of water droplets impimjimj on a wedge in a two-dimensional 
suj)er sonic flow field with an attached shock wave. The closed- 
form solution yields analytical expressions for the equation of 
the droplet trajectory^ the local rate of impingement and the im- 
pingement velocity at any point on the wedge surface, and the 
total rate of impingement. The analytical expressions are 
utilized in the determination of the impingement of water drop- 
lets on the forward surfaces of symmetrical double-wedge airfoils 
in supersonic flow fields with attached shock waves. 

For a wedge, the results provide information on the effects of 
the drojdet size, the free-stream Mach number j the semiapex 
angle, and the pressure altitude. For the double-wed ge airfoil, 
additional calculations jwovide in formation on the effect of air- 
foil thickness ratio, chord length, and angle of attack. 

The results for the symmetrical double-wedge airfoil are also 
correlated in terms of the total collection efficiency as a function 
of a relative modified inertia parameter. The results are pre- 
sented for the following ranges of variables: free-stream static 
tem 2 )erature, 420 ^ to ^6*0° U; droplet diameter, 2 to 100 microns; 
free-stream Mach number, 1.1 to 2.0; semiapex angle for the 
wedge, l.lJf to 7.97^, and corresponding double-wedge-airfoil 
thickness-to-chord ratio, 0.02 to 0 . 14 ; pressure altitude, sea 
level to 80,000 feet; and chord length, 1 to 20 feet. 

INTRODUCTION 

The prol)]em of ice prevention on aircraft flying at sub- 
sonic speeds up to critical Mach numliers has lieen a subject 
of considerable study and research by the NACA. The recent 
advent of aircraft flying at transonic and supersonic speeds 
has required an extension of these icing studies. That an 
icing problem exists in the transonic and sujiersonic speed 
range is verified in reference 1, which shows by an analytical 
investigation with experimental confirmation that diamond 
or symmetrical double-wedge airfoils are subject to possible 
icing at flight Mach numbers as high as 1.4. A similar re- 
sult is expected for other airfoil shapes being considered for 
use at transonic and supersonic flight speeds. 

In conducting research on the ])roblem of ice prevention 
on aircraft and missiles, regardless of tlie magnitude of the 
flight speed, it is essential that the impingement of cloud 
droplets on airfoils and other aerodynamic bodies be deter- 
mined either by theoretical calculations or experiment. The 


impingement variables that must be determined are the total 
water catch, the extent of impingement, and the rate of im- 
pingement per unit area of body surface. These variables 
can be determined analytically from calculations of the cloud- 
droplet trajectories obtained for the various aerochmamic 
bodies. Investigators have reported the results of studies of 
clouddroplet trajectories about right-circular cylinders (refs. 
2 to 5) and about airfoils (refs. 6 to 9) immersed in an in- 
compressible fluid. An evaluation of the effect of compres- 
sibility on the droplet trajectories about cylinders and airfoils 
up to the critical flight Mach number is presented in refer- 
ence 10. 

At present, little information exists on the impingement 
of droplets on aerodynamic bodies in a supersonic air stream. 
The concentration of past effort on problems of impingement 
on aii*foils at subsonic flight speeds and the present lack of 
convenient and rapid means for obtaining the rotational 
flow fields about airfoils at supersonic speeds are possible 
explanations for the scarcity of trajectory calculations for 
the supersonic region. An initial contribution to the solu- 
tion of the over-all problem of imi)ingement of water droplets 
on aerodynamic bodies at supersonic speeds is given in 
reference 11, which presents an analysis of the water- 
interception characteristics of a wedge in a supersonic flow field . 

The present report extends the analysis of reference 11 
and further presents an extensive study of the impingement 
of water droplets on two-dimensional wedges and double- 
wedge airfoils for supersonic flight speeds that result in 
attached shock waves and constant velocity fields behind the 
shock waves. For the wedge angles and double-wedge- 
airfoil thickness ratios considered herein, the shock-wave- 
attachment Mach number varies from a value slightly 
greater than 1 to about 1.4. The method em])loyed is based 
on an analytical solution of the equations of motion by means 
of a closed-form integration. The closed-form solution 
yields anal 3 Tical expressions for the equation of the trajec- 
tories, the local impingement efficiencies, the velocity' at an^^ 
point on the trajectories, and the total rate of impingement. 
This solution is made possible by the use of an empirical 
relation for the drag coefficients for s])heres that gives a good 
ap])roximation of the experimental drag coefficients. 

The results of calculations for the rate, the extent, and the 
distribution of impingement of water droplets on wedges 
and sAunmetrical double-wedge airfoils are presented herein. 


Supersedes NACA TN 2971, “Impingement of Water Droplets on Wedges and Diamond Airfoils at Supersonic Speeds,” by John S. Serafini, 1953. 
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The ranges of variai)les inelnded for the wedge are: free- 
st ream static temi)eratiire, 420°, 440°, and 400° K; dro])let 
diameter, 2 to 100 microns; free-stream Mach number, 1.1 
to 2.0; tangent of the semiapex angle, 0.12 to 0.14; and 
])ressure altitude, sea level, 15,000, and 20,000 iVct. The 
ranges of variables for the doubk'-wedge airfoil are the same 
as those' for the wedge, and the additional variables for the 
symmetrical double-wedge airfoil range from 0.02 to 0.14 for 
the thickness ratio and from 1 to 20 feet for the chord lengU\. 

The work n'jiorted herein was jierformed at the' XACA 
Lc'wis laboratory in the sjii’ing and summer of 1952. 

ANALYSIS 

STATEMENT OF PROBLEM 

The solution of the problem of impingement of water 
droplets on a two-dimensional wedge at supersonic speeds 
with an attached shock wave is not as difficult as that for 
the impingement on various airfoils at low subsonic speeds. 
For the wedge at supei'sonic speeds with an attached shock 
wave, the air velocity everywhere ahead of the shock wave 
is constant and ecpial to the free-stream air velocity Vj 
(fig. 1). The air velocity behind the shock wave is also 
everywhere constant and i)aj’allel to the wedge surface. All 
the droplets have tlie same initial velocity (that of the free- 
stream air velocity), and their trajectories are exactl}' co- 
incident with the air streamlines upstream of the shock wave. 
All water droplets of a given size are subjected to identical 


air-velocity fields, whicii in turn ])roduce identical force 
systems downstream of the shock wave, irresi)ective of the 
])oint along the shock wave where the droplets cross the 
wave. It follows, therefore, that, for droplets of a given 
size, all the trajectories in a givT'ii ])rohlem are identical with 
respect to the point where the droi)let crosses the shock wave. 
By ado]3ting a frame of reference that moves at the con- 
stant velocity of the air ^2 downstream of the shock wave, 
the prol)lem of the drojilet motion is reduceil to the still-air 
])rohleni, defined as the determination of the motion of a 
dro])let that, having an initial velocity, is j)rojected into 
([uiescent air. Hence, relative to the moving frame of refer- 
ence, the initial velocity of the droplet upon crossing the 
shock wave is ('(pial to the vectorial difrerence of the free- 
stream air velocity Vi and the air velocity 1^2 downstream of 
the shock wave. Adoption of the frame of refen'iice moving 
with a constant velocity reduces the problem from the solu- 
tion of two simultaneous nonlinear second-order differential 
e([uations in the fixed coordinate system to the solution of a 
single nonlinear second-order differential equation in the 
moving cooi-dinate system. At any given instant, the drop- 
let dis])lacenient relative to the point of intersection with the 
shock wave in the fixed frame of reference is obtained by 
adding vectorially the droplet displacement within the mov- 
ing frame of reference to the displacement of the moving 
frame of reference for the same increment of time. 


At /^O At /^/ 
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This general nietliod was used in reference H, where the 
one second-order difi’erential equation representing the drop- 
let motion relative to the air velocity l)ehind the shock wave 
was integi’ated graphicallj^. IToweve]*, it is possi))le to 
ol)tain a com])letely analytical solution by means of a closcal- 
form integration without resorting to the use of numerical 
integrations or analog computing equipment, if an empirical 
relation is assumed for the drag coefficient as a function of the 
Reynolds number of the droplet relative to the air. It will 
be shown that this closed-form integration of the still -air 
problem wlien applied to the wedge in supersonic flow with 
attached shock wave yields the equations for the trajectories 
of the water droplets and the di'oplet velocities at any point 
on the trajectories and makes available relations for the 
rates of total water impingement and the local rates of water 
impingement along the wedge surface. Fuz-thermore, it is 
shovTi that these equations can also be readily applied to the 
determination of the droplet impingement on a double-wedge 
aii’foil in supei’sonic flow vvdth attached shock wave. 

Three of the usual assumptions made in the previous in- 
vestigations on impingement at subsonic s])eeds and also 
required for this investigation az-e (1) the watez* dz-oplets az*e 
always spherical and do not change in size, (2) the force of 
gi-avity on the droplet may be neglected in coznpaz'ison with 
the dz*ag forces, azzd (3) the drag of the air ozi the dz-oplet is 
that of a viscous incompressible fluid. Hei-e it is additionally 
assumed that (4) the two-dimensional supersonic flow field 
about the wedge is fz'ictionless except within the infinitesi- 
mally thin attached shock wave, (5) no condensation shock 
occurs and no change in phase occurs as the water dro])lets 
tz-avez*se the oblique shock wave, azid ((>) the uzibalance of 
the forces on the water di-oplet from the instant it enters the 
shock wave until it emerges fi*om the shock wave can be neg- 
lected izi the calculation of the trajectozles. 


EQUATION OF DROPLET MOTION IN MOVINC REFERENCE FRAME 

The velocity of the dz-oplet in the znoving fz-ame of z-efez*- 
ence is 

U=\Vh~V2\ ( 1 ) 

where V ^ is the droplet velocity with respect to the fixed 
fi-ame of reference, and I "2 is the air velocit}" downsti-eam of 
the attached shock wave also with respect to the fixed 
fi-ame of reference (fig. 1). In the fi*ame of z-eference moving 
with the velocity the statement of Newton’s law of mo- 
tion for the water droplet becomes 

- Cz, I ( 2 ) 


from which 


d(7 3 p2 ^ 


(3) 


(A complete list of s^unbols is given in appendix A.) 

Equation (3) is the differential equation of motion of a 
di-oplet pi’ojected with an initial velocity into a region of 
quiescent air (the so-called still-air problem). The shock 
wave is considered to be a surface of discontinuity from which 


the di*oplets emei-ge with a velocity TV I^i this case the 
initial velocity of the droplet is 


Cr-jv^-V^I (4) 

which is the magnitude of the vector difference of the air- 
velocity vectoi-s upstream and downstream of the attached 
shock wave. As can be shown fi-ozn considei-ation of the 
continuity equation and the equation for conservation of 
momentum aci-oss the oblique shock wave, the velocity vec- 
tor Ui is normal to the shock wave. At any subsequent 
instant of time, tlie relative di-oplet velocity vector U retains 
the same angular oi-ientation to the shock wave and changes 
only in magnitude. 

In refei-ence 11 the solution of equation (3) is obtained by 
numei’ical integration. The result obtained in this manner 
makes it necessaz-y to use a gz-a])hical pz-oceduz*e izi determizi- 
ing the tz-ajectories and the local rates of iznpizzgeznent. 
Howevez-, an analytical solution of equation (2), which elimi- 
nates the gz-aphical pz-oceduz-e, cazi be obtaizzed if the experi- 
znental values of the di-ag coefficient Cd az-e expressed in a 
function involving the Reynolds zzunzbor The z*elation 

is 

= (1 + (5) 


where e azid m are the eznpiz’ical cozistants. This eznpiz’ical 
relation is a valid appz-oxiznatiozi in the z-azzge of Reynolds 
ziuznbez-s to which cloud dz-oplets are subjected in trajectory 
calculatiozis. Substitution of the expressiozi for Cd (oq. (5)T 
in equation (3) results in the expressiozi 


(It (If 2 pu, 



( 6 ) 


where the local relative RoAuiolds ziuznber Ber = 2p2UalpL2- 
The displaceznezit of the water dz-oplet in the znoving fz-ame 
of 1 -efez‘once x is measuz-ed fz-om the air stz-eamline that intez-- 
sects the shock wave at the poizit where the water dro])let 
enters the air-flow field dowzistreazn of the shock wave. 
The closed-foi’in integi-ation of the diffoz-ential equation (6) 
is jiresezited in appezidix B. The use of 2/3 for the expoziezit 
m and 0.158 for the value of the eznpirical constant e in equa- 
tion (5) yields an einpii-ical cuz-ve for the dz-ag coefficient as 
a functiozi of the local Reyziolds nuznbers that approximates 
very well the variatiozi of the experimental values of the 
dz-ag coefficient in the z-ange of Reynolds mzmbez-s fz-om 0.5 
to 500. The value of 2/3 for the expoziezit m also facilitates 
the closed-foz-zn integz-ation of the differential eqization of 
motion. In figuz-e 2 a graph of the empiz-ical relation is 
])z-esented, along with the drag-coefficient data of referezices 
4 and 12. 

The rcsizlts of the integz-ation az-e given by tlic following 
equations: 


3 P2 


2 L 


( 7 ) 


4 


REPORT n.')9 NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


1 ) 6 ^- 1 ]- 3 /^ ( 8 ) 

and 

wIkto 

Tlio intonmaliate slops of into^ralion aro given in a])])ondix B. 

Equations (7) and (8) give, ros])ootivoly, the dis])laeement 
and the veloeity of the dro])let at any instant in the moving 
frame of referenee. The displacement of the dro])let with 
res])ect to the ])oint wliere it crosses the shock wave can he 
obtained i)y a vectorial addition of the displacement x and 
the dis])lacement of the moving reference frame in the cor- 
res])onding time interval. The dro])let velocity relative 

to tlie fixed frame of reference must also be obtained by the 
vectoi'ial addition of [' (ecp (8)) and Ecpiation (9) 

gives the maximum value of x obtained as the time of travel 
in the air^llow field downstream of the shock wave 
ai)])roaches infinity. The significance of this (juantity will 
be discussed in s\d)se(pient sections. 


REL.\TIONS RE(iUIREI) FOR .\FFLICATION OF ('LOSEI)-FORM SOLUTION 
TO OBTAIN DROFLET MOTION AND IMFINGEMENT IN FIXED 
REFERENCE FRAME 

Impingement on wedges. For a ])i-oblem of given aero- 
dynamic conditions, the trajectories of all the water di‘ 0 ])lets 
of a given size are identical when the points where the droj)let 
t rajectoi’ies intersect the shock wave are su])eriniposed. 
This unique characteristic of the water-drojilet trajectories 
about a wedge in sujiersonic flow with an attached oblique 
shock wave is the result of two constant velocity fields, one 
u])stream and one downstream of the shock wave. There- 
foi-e, only one set of equations for a single trajectory is neces- 
sary to calculate the im])ingement ])aranieters for a specified 
problem, including a given droj)let size. The values of the 
initial relative velocity Ut, the initial Reynolds number 
Rer,i, and the density ratio p^^/po are needed for substitution 
in the closed-form solution of the ecpiations of droplet motion. 
These values can be obtained from information available in 
reference 13 and fi'om the use of simple algebraic and ti*igo- 
nometric relations for given values of the free-stream static 
temperature h, the dro])let diameter d, the free-stream Mach 
number d/i, the angle of surface inclination to the free-stream 
direction a, and the free-stream static, pressure /q. These 
relations result in the following exjiressions for the initial 
relative velocity and the initial relative Reynolds number: 



Droplet Reynolds number, Re^ 


Fin V RE 2. — Comparison of empirical relation for drag coefficient as function of droplet Ueyuolds number with experimental values. 


Colder 
24 ~ 


l + 0.158/eer^/« 
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and 


[/,= |C/,|=FiVl+f2"-2Q cos (T*= V’.co (10) 

( 11 ) 

\ git Pi 


A coiiveilieiit form of the solution for the impingement on 
a wedge or the front half of a double-wedge airfoil is obtained 
if S is defined as the distance to the point of impingement 
measured from the leading edge for a water droplet that 
enters the flow field behind the shock wave at a distance f 
above the leading edge (fig. 1). The unique relation between 
S and f in a given problem for droplets of the same size is 
quite readily determined by considering the displacement of 
the water droplets as the vectorial sum of the displacement 
of the water droplet relative to the moving frame of reference 
and the displacement of the moving reference frame relative 
to a fixed frame of reference (referred to wedge). Since the 
moving reference frame has a velocity equal to the air velocity 
T" 2 , which is constant in magnitude and parallel to the wedge 
surface, only the droplet travel in the moving reference frame 
includes the component of droplet travel representing the 
approach of the water droplet to the wedge surface. For a 
water droplet starting from j^oint A and impinging on the 
wedge surface at point D (fig. 1), the displacement of the 
moving reference frame (=V 2 t, where t is zero at p oint A) is 
given by the displacement vector AB equal to CD, and the 
droplet motion in the moving frame of reference is given by 
the displacement vector BD equal to AC. Therefore, rela- 
tive to the starting point at A (fig. 1), the displacement of 
the water droplet to the point of i mpin gement at D is 
obviousl}^ equal to AB + BD or to AC + CD. From figure 1, 
the displacement of the water droplet at the point of impinge- 
ment D, measured from the leading edge at E, is given by 
adding the vector EA to the displacement vector from the 
starting point A; and the displacement of the droplet at D 
referred to the leading edge is 


S=lEDl = |EA + AC + CDl = lEC + CD| = |EC| + |CD| = ^+r 

( 12 ) 


where ^ is the magnitude of the displacement vector EC, 
and i=V 2 t) is the magnitude of the displacement vector 
CD (the displacement of the moving frame of reference). 

The values of f and ^ are obtained in terms of x, the distance 
of travel in the moving reference frame, from sinij)le trigo- 
nometric identities involving the various angles shown in 
figure 1: 

sin B tan (13) 

where x is given by equation (7); r by 

-1 (13a) 

CajJ 

and 

^=x sec (p-^a) (14) 



Substitution of equation (14) into equation (12) for the sur- 
face distance to the point of impingement yields the following 
equation: 


2 

S = x sec -h V 2 T (1 5) 

0/i2 

Since a: is a function of r in equations (13) and (15), the 
exi)ressions for f and S, respectively, are functions of r. 
However, since r cannot be eliminated from equations (13) 
and (15), S cannot be obtained explicitly as a function of f. 
Nevertheless, the curves of f, the initial displacement of the 
water-droplet trajectory from the leading edge normal to 
the free-stream direction, against *S, the distance to the point 
of impingement of the stated water-droplet trajectory, can 
be obtained by substitution of the same set of values for r 
in the expressions for f and S. 

An analytical expression for the local impingement 
efficiency can be obtained from the preceding expressions 
for f and S. The local impingement efficiency jS is defined 
by the expression 

(Jt 


where Af is the difference in the initial displacements of two 
water droplets having very nearly equal initial displacements, 
and AS is the small increment of wedge surface between the 
points of impingement of the two water droplets. Differen- 
tiating f and S with respect to r and performing the division 
indicated by equation (16) yield the following expression 
for i8: 


where 




1 


ni Ui a 


Pw 


(17) 


rii = sin 6 tan (r + cr) 


(18a) 


712= sec 


(18b) 


773 — 1^2 


PwCl^ T- t 
^ \ 2 - 
0/X2 T 


(18c) 


Since ^ and S' both are functions of r, the local impingement 
efficiency at any point on the wedge surfaee is determined 
by using the same value of r in equations (17) and (15). The 
value of jS that exists as the point of impingement of the 
water droplet on the wedge surface approaches the leading 
edge as a limit (S-^0) is given by the following: 


/3o=lim (13) 


1 

^ ^ J_ 

Til 77 1 Pu, IJ i 


sin a 


(19) 


The magnitude and direction of the droplet velocity at the 
l)oint of impingement Va,im (relative to the fixed frame of 
reference) can also be easih^ obtained at any point on the 
wedge surface as 

Va.i,n= F. (-^y-t-2Qo, ^ cos (v + a) (20) 

/c,„ = «r — tr'„=cr— sin"* sin (i"f<r)1 (21) 

where k is the angle between the free-stream velocity vector 


G 
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Shock wove Expansion zone 





(c) 


and (he droplet-velocity vector ^^id a' is the angle 
between the droplet-velocity vector and (he air-velocity 
vector VV In equations (20) and (21) for a given trajectory, 
Ti, 0, CO, Ui, V, and a are constants. Therefore, for a given 
trajectoiy, Va,im ‘^nd Kim arc functions only of U (eq. (8)), 
which is in turn a function of r, the dimensionless time 
variable. 

Impingement on symmetrical double-wedge airfoils. — 
The impingement on a double-wedge airfoil may he oldained 
from the solution to the problem of impingement on a wedge 
as presented heretofore. In this report, the double-wedge 
airfoil considered is symmetrical, the maximum thickness 
occurring at 50 percent of chord (fig. 3). At zero angle of 
attack, the impingement on a double-wedge airfoil will be 
limited to the region from the leading edge to the shoulder 
at 50 i)ercent of chord. The solution for impingement on a 
wedge surface having a given semiapex angle a (the angle of 
inclination of cither wedge surface to the free-stream direc- 
tion) can also be used as the solution for a doul)le-wedge air- 
foil where the thickness ratio is equal to tan cr, the tangent of 
the semiapex angle, and where the droj)let size and other 
j)aranieters of the ])robleni are the same as for the wedge. 
Therefore, the values of the local impingement efficiencies 
(3 and j3o at any given point on the surface will be identical 
for both the wedge and the double-wedge airfoil at zero angle 
of attack under the aforementioned similai*ity of conditions. 

The solution for the impingement on the double-wedge 
airfoil at angle of attack can also be obtained from the 
solution for impingement on a wedge as for the case of the 
double-wedge airfoil at zero angle of attack. AVlien the 
symmetrical double-wedge airfoil is at angle of attack a, 
the angle of inclination of its forward u])per surface to the 
free-stream direction is equal to a— a and that of the forward 
lower surface is equal to a-\-a. Therefore, the solution to 
the im])ingemcnt on the u])])er and lower surfaces of the 
double-wedge airfoil is obtained from tlm solutions for 
im])ingement on wedges having the redefined semia])ex angles 
of a— a and (r + a, resi)ectively, where the dro])let size and 
other parameters of the problem are ke])t the same. For 
the double-wedge airfoils at angles of attack having tangents 
equal to or greater than the thickness ratio, the water 
droplets will not impinge on the up])er surface. At angles 
of attack having tangents greater than the thickness ratio, 
some water droplets may iin])inge on the lower surface 
beyond 50 percent of chord. These three conditions are il- 
lustrated schematically in figure 3. For of<tan"^(T/c), the 
im])ingement occurs on surfaces AC and AB;for «= tan” ^(7"/c), 
impingement occurs only on surface AC; for a)>tan” ^(77c), 
ini])ingement occurs on lower surface AC and may occur 
on lower surface CD. However, the condition where 
a^tan” '(T/c) is not considered herein, since the solution 
presented in this ]*e])ort is not valid for the determination 

(a) Angle of attack ct<C\siU~^{Tlc) . 

(b) Angle of attack a = tan~qT'/c). 

(c) Angle of attack a^taiw^T/r) . 

Figure 3. — Schematic diagram of symmetrical double-wedge airfoil at 

angle of attack in supersonic flow with attached shock wave. 
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of trajectories of droplets iinpiiigiiig; on the surface beyond 
the shoulder or 50 percent of chord of the double-wed^e 
airfoil (surfaces BD or CD), where a portion of the trajectories 
is within the expansion zone emanating from the shoulder. 

RESULTS AND DISCUSSION 

From the equations presented in the previous section and 
in appendix B, the impingement of water droplets on a wedge 
in a supersonic flow field with an attached shock wave can 
he calculated over a range of free-stream conditions, wedge 
angles, and droplet sizes. As has been indicated previously, 
the impingement characteristics of a double-wedge airfoil 
(at zero angle of attack and also for small values of angle ol 
attack) can readily be determined from the impingement on 
wedges for similar conditions. The results for the wedge 
and the double-wedge airfoil are pi-esented and discussed 
se])arately. A comparison of the total collection efficiency 
and the water impingement rate at zero angle of attack for a 
double-wedge airfoil with those for an NACA 0006-64 
airfoil is presented in appendix C. 

WKDGE 

Local impingement efficiency. -The rate of water impinge- 
ment on a local area of wedge or airfoil surface is proportional 
to a dimensionless term jS, the local impingement efficiency. 
The local rate of water impingement in ])ounds ])er hour 
])er square foot is 

where j8 is defined in equation (16). The local im])ingenient 
efficiency, when given as a function of tlie surface distance 
of the wedge, allows the determination of the local rate of 
impingement of water droplets at any ])oint on the surface, 
the total impingement of water droplets on the entire sur- 
face or any given portion of the surface, and the extent of 


impingement on the surface. The local impingement 
efficiency is related to a point at a given distance S on 
the wedge surface in equations (15), (17), and (19) by the 
dimensionless time variable r, which is common to all three 
expressions. 44ie variation of 0 with S is pi*esented in 
figure 4 for an extensive range of free-stream conditions, 
semiapex angles, and droplet sizes. The value of ^ at the 
leading edge (aS-^O) is the sine of the semiapex angle (sin a); 
and as aS increases, the value of 13 decreases rapidly and ap- 
proaches the value of zero asymptotically as S api)i*oaches 
infinity. However, it is to be noted that negligibl 3 " small 
values of ^ (/3~1 percent of for the wedge) are attained 
at large but finite values of S. 

The curves of as a function of S presented in figure 4 
are those of an idealized situation. The assumed two- 
dimensional supersonic flow field about the wedge does not 
account for a stagnation point that must exist at the leading 
edge of the wedge, regardless of the sharpness of the leading 
edge. In addition, the leading edges of wedges and double- 
wedge airfoils might be considered to be somewhat rounded 
when subjected to considerable magnification. Therefore, 
it is reasonable to assume that very near the leading edge 
(S~^0) the value of p would actuall^^ be greater than the 
calculated value of at the given distance S. However, 
this should have a negligible effect on the rest of the incurve 
and also on the total impingement on the wedge, since the 
effect of a stagnation jioint Avould be limited to a veiy 
small region about the leading edge. 

For S approaching veiy large values, the calculated values 
of probably difl’er somewhat from actual values obtained 
in flight, because the anah tical solution of the present report 
does not consider the existence of the boundary la^TT on the 
wedge surface. Since the boundaiy-lay^er thickness increases 
with the surface distance along the wedge, droplets that ini- 



(a) Effect of pressure altitude. Droi)lct diameter, 20 microns; free-stream Mach number, 1.4; tangent of semiapex angle, 0.02. 
Figure 4. — Local impingement efficiency on wedge as function of distance along surface. Free-stream static temperature, 440° R. 
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pingo at large distances from the leading edge actually would 
have traveled through the boundarv layer for some non- 
negligihle interval of time. However, only a very small 
fraction of the total water droplets of a given size impinge 
under this condition, and for large values of S the values of 
/3 are negligibly small. For example (fig. 4 (a)), such would 
he the case for values of greater than 8 or 9 feet, where 
/3<0.0002 as com])ared with ft, ^0.020. 

A preliminary survey disclosed a negligible effect of the 
free-stream static temperature on the local imi)ingement 
efheieney as a function of the surface distance along the 
wedge {13 against S). Values of (3 and corresponding values 
of S were calculated for free-stream static temperatures of 
420°, 440°, and 460° 11, droplet diameter of 20 microns, free- 
stream Mach number of 1.3, tangent of the semiapex angle 
of 0.06, and ])i*essure altitude of 15,000 feet. The values of 
p for the free-stream static temperatures of 420° and 460° K 
are within 1 percent of the values of at 440° K. Since these 
calculations show that curves of ^ against S for the three 
values of free-stream static temperature form ])ractically a 
single curve when ])lotted to the usual scales, no figure is 
presented to illustrate the effect of free-stream static tem- 
perature on imiiingement. Furthermore, the results included 
herein, which are calculated for a free-stream static temper- 
ature of 440° K, may he used in the range of temperature 
from 420° to 460° K or possibly an even greater range. 

The effect of the free-stream static ])ressure on the local 
im])ingement efficiency as a function of the surface distance 
along the wedge is ])resented in figure 4 (a) for ])ressure 
altitudes of sea level, 15,000, and 30,000 feet. Increasing 
the pressure altitude (decreasing the free-stream static ])i*es- 
sure) increases the values of /3 at any distance S. For ex- 
ample, at 8=1.5 at sea level, )3 is 0.0052; and at 30,000 feet, 
/3 is 0.0071. Since po (the value of at 8=0) is ecpial to 
sin O', where a is the semiaiiex angle, the curves for the three 
pressure altitudes have the same maximum value of Also, 
for the various ])ressure altitudes, the extent of impingement 
along the wedge surface is essentially the same. 

The effect of the semiapex angle o- of the wedge on the 
local impingement efficiency as a function of the surface dis- 
tance along the wedge is ])resented in figure 4 (h) for values 
of tan a from 0.02 to 0.10. Since the values of ft, are equal 
to sin O', increasing the semiapex angle of the wedge results 
in an increase of ft,. The surface extent of perce])tihle im- 
jiingement (as characterized i)y /^^^^O.Olft,) does not vary as 
the wedge thickness is increased. 

The effect of free-stream Mach number on the 13 curve is 
presented in figure 4 (c) for free-stream Mach numbers of 
1.2, 1.3, 1.4, 1.5, and 2.0. For the wedge semiapex angle 
presented in the figuiv (tan o- = 0.04), the value of J/i=1.2 is 
close to the shock-wave-attachment Mach nuniher. The 
shock-wave-attachment Mach nurnher is a function of the 
wedge semiapex angle and is defined as that Mach nurnher 
below which the shock wave is detached from the wedge. 
An increase in the free-stream Mach number results in 
an increased surface extent of jierceptible impingement and 
also in an increased value of at any given distance 8 (exce])t 



(b) Kffc'ft of wedge semiapex angle. Droplet diameter, 20 microns; 
free-stream Mach number, 1.3; pressure altitude, 15,000 feet. 

Figure 4. — Continued. Local impingement elliciency on wedge as 
function of distance along surface. Free-stream static temperature, 
440° U. 

at 8=0, whore |3 = 0o = sin cr and at 8-^oo, where j3->0). 
This increase in the surface extent of ])ercej)til)le imjiinge- 
ment is shown in figure 4 (c), in which, for free-stream Mach 
numhers of 1.2, 1.3, 1.4, 1.5, and 2.0, the surface extents 
of ])erce])til)le imjiingement on the wedge (where 0.01 ft)) 
are 5.35, 6.05, 6.65, 7.20, and 9.4 feet, respectively. 

The effect of the droplet size on ^ is presented in figure 4 (d) 
for droplet diameters of 10, 20, 30, 40, 50, and 100 
microns. The surface extent of impingement and the values 
of P at any given distance 8 are considerably increased as the 
droplet size is increased. For example, for the semiapex 
angle presented in the figure (tan (r=0.06), at 8=3 feet the 
values of 0 are 0.0000, 0.0033, 0.0106, 0.0182, 0.0244, and 
0.0417 for values of droplet diameter of 10, 20, 30, 40, 50, 
and 100 microns, respect ivel}L The surface extent of per- 
ceptil)le impingement has values of 1.5, 5.7, 11.2, 18.1, and 
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(c) Effect of free-stream Mach number. Droplet diameter, 20 microns; tangent of semiapex angle, 0.04; pressure altitude, 15,000 feet. 

(d) Effect of droplet size. Free-stream Mach number, 1.3; tangent of semiapex angle, 0.06; pressure altitude, 15,000 feet. 

Figure 4. — Concluded. Local imj)ingement efficiency on wedge as function of distance along surface. Free-stream static temperature, 440° R. 
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28.0 feet for droplet diameters of 10, 20, 30, 40, and 50 
mierons, respectively. As shown by the preceding dis- 
cussion and also hy a comparison of the values of (3 as a 
function of S (fig. 4), varying the droplet diameter from 20 
to 30 microns or from 30 to 40 microns is of the same order 
of magnitude in its effect on the against S curve as varying 
the ])ressure altitude from sea level to 30,000 feet or varying 
the free-stream Mach number from 1.2 to 1.5. 

Total impingement. The effects of the free-stream Mach 
iiumher, the semiapex angle of the wedge, the pressure alti- 
tude, and the dfoplet size on the total impingement on a 
wedge surface of infinite extent are given in figure 5. The 
total im])ingement is re])resented by which is the f of 
tlie droplet having its trajectory tangent to the wedge surface 
(theoretically the tangent trajectory touches the wedge 
surface at a point aV— > oo ). This can he obtained from the 
integration 

^ 00 

= fids 

or more directly from ecpiation (13) of the analytical solution, 
as 

= sin 0 tan (v+a) 

where x,n is given hy eciuation (9) and is also defined in api)en- 
dix A. The value for is obtained from the expression for 
-c (eq. (7)) by allowing r^oo. Since the droplet that enters 
the flow field downstream of the shock wave at a distance 
^ ecpial to (fi^- 1) theoretically impinges on the wedge 
surface only as the surface distance S ap])roaclies infinity, 
only di’oplets having values of ^ equal to or less than ^rn ^vill 
inij)inge on the wedge surface of infinite extent. The rate 
of total water catch on one wedge surface in terms of is 
given as 

U;, = 0.329()f,A pc 

where II is expressed in ])ounds per hour ])ei* foot of span, 
in miles ])er hour, and iv in grams ])er cubic meter. 
Therefore, is directly proportional to the rate of total 
water catch on the entire wedge surface and is the rate of 
total water catch on one wedge surface ])er unit of span, 
free-stream velocity, and liquid-water content (in appropri- 
ate units). 

The variation of with the tangent of the wedge semia])ex 
angle tan a is shown in figure 5 (a) for free-stream Mach 
numbers of 1.2, 1.3, 1.4, and 2.0. As ex])ected from the 
curves of ^ as a function of S, the value of increases as 
tan a increases. However, tlu' rate of incj*ease in with 
res])ect to tan a decreases with an increase in tan a. As can 
be seen from figure 5 (a), increasing the free-stream Mach 
numb('r increases the value of for a constant value of tan a. 

The variation of f,,, with ])ressure altitude is ])resented in 
figure 5 (b) for two Mach mimbers. In the range of ])ressure 
altitude from sea level to 30,000 feet, the increase of with 

(a) Effect of tangent of seiniapex angle. Droi)let diameter, 20 microns; 
])ressnre altitude, 15,000 feet. 

(h) Effect of pressure altitude. Droplet diameter, 20 microns. 

(c) Effect of droplet size. Pressure altitude, 15,000 feet. 

Eigure 5. — Total impingement rate ir,„ for wedge of infinite extent. 

Eree-stream static temperature, 440° 11; ir,„ = 0.3296^„iEii«^- 
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an increase in pressure altitude is ap])roxiinately linear. 
Tile variation of witli tlie droplet diameter d in microns is 
shown in figure 5 (c). In the range of droplet diameter from 
10 to 100 microns, ^rn as a function of d results in a curve 
that is very nearly a straight line when plotted on logarithmic 
paper. This linearity permits an accurate interpolation of 
when calculations are made for a few droplet diameters 
for a given value of free-stream Mach number, wedge semi- 
apex angle, and pressure altitude. 

Droplet velocities at impingement.— The variation of 
Irf.rw/li (ratio of droplet impingement velocity to free- 
stream velocity) with the surface distance along the wedge 
is presented in figure 6 for three cases. These three cases are 
representative of the results when the droplet diameter d is 
20 microns and the pressure altitude is 15,000 or »10,000 feet. 
The curves of i as a function of S have characteristics 

similar to the curves presenting as a function of *S. At 
S={), obviously, all the curves have equal to unity; 

and, as S is increased, the value of the velocity ratio rapidly 
decreases and asymptotically approaches V 2 IV 1 , the ratio of 
the air velocity downstream of the shock wave to the air 
velocity upstream of the shock wave. The solid curve in 
figure 6 illustrates a typical situation for which Vyi'i is very 
nearly unity (tan o-=0.02 and Mi = 2A)). The two lower 
curves (tan o-=0.06 and 0.10 at A/i=1.3) are typical for cases 
where a stronger shock wave produced by a larger semiapex 
angle results in decreased values of for large values 

ofS. 

SYMMETRICAL DOUBLE-WEDGE AIRFOIL 

As shown in the ANALYSIS, the local impingement 
efficiency ^ at any point on the forward surfaces of a double- 
wedge airfoil (surfaces AB and AC in fig. 3) can be obtained 
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Figure G. — Variation of ratio of impingement velocity of droplet to 
free-stream velocity with distance along surface of wedge. Free- 
stream static temperature, 440° R; droplet diameter, 20 microns. 

directly from the results for the ^ against S curves for wedges. 
The local impingement efficiency ^ may be obtained from 
figure 4 or equation (17). 

In general, the results for the impingement on a sym- 
metrical double-wedge airfoil are presented in this report in 
terms of the total collection efficiency as a function of the 
scale parameter in an attempt to conform with the 
existing literature on the impingement characteristics of 
airfoils. In the notation of the present report, the total 
collection efficiency Em as stated in references 7 and 8 
is defined as 

= + ( 22 ) 



(a) Effect of pressure altitude and chord length. Droplet diameter, 20 microns; free-stream Mach number, 1.4; airfoil thickness ratio, 0.02. 

Figure 7. — Total collection efficiency of symmetrical double-wedge airfoils as function of scale parameter. Free-stream static temperature, 

440° R; angle of attack, 0°. 


292971— .54 3 



12 


REPORT llo9 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


wIkmt T is the inaxinuim thickness of the symmetrical 
(loiihle-wedge airfoil, and and f/ are the absolute values of 
tlie initial displacements from the leading edge (in a direction 
normal to free-stream direction) of the droplet trajectories 
that impinge at the sho\ilder of the upper and lower sur- 
faces, respectively, of the double-wedge airfoil. For the 
symnu'trical double-wedge airfoil at zero angle of attack, 
and if,j will be equal; at an angle of attack, the tangent of 
which is c(pial to the thickness ratio, the value of is equal 
to zero. 

scale parameter xp is calculated for the double-wedge 
airfoil at supersonic speeds as for other airfoils at subsonic 
speeds. It is defined as 


where c is the chord length of the double-wedge airfoil. The 
results presented for impingement on double-wedge airfoils 
use essentially the same parameters as used for impinge- 
ment on wedges — droplet size, free-stream Mach number, 
double-wedge-airfoil thickness ratio, and pressure altitude. 
In addition, angle of attack and chord length are specified. 

Total collection efficiency at zero angle of attack. -The 
variation of the total collection efficiency Ern ^vifh the scale 
parameter \p is presented in figure 7 for a symmetrical 
double-wedge airfoil at free-stnTim static temperature of 
440° K and zero angle of attack. The effect of pressure 
altitude on the variation of E^ with \p for a symmetrical 
double-wedge airfoil of 0.02 thickness ratio is given in figure 
7 (a) for dro])let diameter of 20 microns and free-stream Mach 
number of 1.4. The ])ressure altitudes presented in the 
figim' are sea level, 15,000, and 30,000 feet. The lines for 
constant values of chord length from 1 to 20 feet are also 
included in the figure. For the double-wedge airfoil sub- 
jected to a constant-velocity supersonic flow field, the total 
collection efficiency E,n increases slightly as the pressure 
altitude is increased when the chord length and the other 
variables are held constant. Considering E„, as a function 
of scale parameter xp where the di’ 0 |)let free-stream Reynolds 
number Hei is held constant yields results similar to those 
for airfoils having rounded leading edges at subsonic s])eeds. 
With Rex held constant, the value of E^ decreases as \l/ in- 
creases in the manner indicated in figure 7 (a). As pre- 
viously mentioned, the total collection efficiency and the 
impingement rate at zero angle of attack for a symmetrical 
double-wedge airfoil are compared in a])pendix C with 
thos(» for an XACA ()()0()-64 airfoil. 

The effect of airfoil thickness ratio on the variation of E,n 
with \p is ])resented in figure 7 (b) for thickness ratios from 
0.02 to 0.14, droplet diameter of 20 microns, free-stream 
M ach number of 1.4, and pressure altitude of 15,000 feet. 
The droplet free-stream Reynolds number is maintained 
at a value of 453 for all the curves. The effect of increasing 
the airfoil thickness ratio is to decrease the total collection 
efficiency. The rate of decrease in Em with an increase in 
airfoil thickness ratio becomes somewhat smaller as the 
thickness ratio increases. 

The effect of free-stream Mach number on the variation 
of total collection efhciency Em with scale parameter xp is 


shown in figure 7 (c) for free-stream Mach numbers of 1.1, 

1.2, 1.3, 1.4, 1.5, and 2.0, droplet diameter of 20 microns, 
airfoil thickness ratio of 0.02, and pressure altitude of 15,000 
feet. The results show that the total collection efficiency 
increases as the free-stream Mach number increases. Ffow- 
ever, the increase in the total collection efficiency from the 
free-stream Alach number of 1.1 to 1.2 is considerably 
greater than the increase in efficiency from a Mach number of 
1.2 to 1.3 and from 1.3 to 1.4, and so forth. Calculations 
for thickness ratios of 0.04 and 0.00 indicate the same trend. 
P"or the 0.02 thickness ratio, the lowest Mach number pre- 
sented in figure 7 (c) (d/=l.l) is quite close to the limiting 
Mach number for shock-wave attachment. Therefore, the 
rate of decrease of total collection efficiency with decreases 
in free-stream Mach number increases as the Mach number 
approaches the shock-wave-attachment Mach number as a 
limit. 

The effects of droplet size and chord length on the varia- 
tion of Em with xp are shown in figures 7 (d) and (e). Figure 
7 (d) presents curves of Em against xp for a symmetrical double- 
wedge airfoil of 0.02 thickness ratio at a free-stream Mach 
number of 1.4 and a pressure altitude of 15,000 feet. The 
curves are for constant values of droplet diameter {d= 
10, 20, 30, 40, and 50 microns) as well as for constant values 
of chord length (c=l, 2, 4, 8, and 20 ft). Increasing the 
dro])let size greatly increases the total collection efficiency. 
For example, at c = S feet, the values of Em are 0.096, 0.310, 
0.495, 0.625, and 0.711 at droplet diameters of 10, 20, 30, 
40, and 50 microns, respectively. The rate of increase in 
the total collection efficiency as the droplet diameter in- 
creases is less for the larger dro])let sizes. This effect can 
also be observed in hgure 7 (e), in which curves of Em as a 
function of xp are presented for a symmetrical double-wedge 
airfoil of 0.06 thickness ratio for droplet diameters of 2, 10, 
20, 30, 40, 50, and 100 microns, free-stream Mach number of 

1.3, and ])ressure altitude of 15,000 feet. 

A coni])arison of figures 7 (b) and (c) (same droplet size 
and pressure altitude) shows that the effect on the total 
collection efficiency Em of an increase in the free-stream 
Mach number from 1.1 to 2.0 is, in general, of the same order 
of magnitude as a decrease from 0.14 to 0.02 thickness ratio. 
For example, for xp of 1790 (c = 8 ft) in figure 7 (b), the values 
of Em decreased from 0.310 to 0.165 for an increase from a 
0.02 to a 0.14 thickness ratio, respectively, a decrease of 
0.145 in the value of E„i. For xp of 1790 (c — 8 ft) in figure 
7 (c), the values of Em increased from 0.230 to 0.400 for an 
increase in the free-stream Mach number from 1.1 to 2.0, 
respectively, an increase of 0.170 in the value of 

Another comparison of figures 7 (c) and (d) (same airfoil 
thickness ratio and pressure altitude) shows that the effect 
on the total collection efficiency Em of an increase in the 
droplet diameter from 10 to 50 microns is much greater than 
an increase in the free-stream Mach number from 1.1 to 2.0. 
For example, for c = 8 feet in figure 7 (d), the value of Em 
increased from 0.095 to 0.710 for an increase in the droplet 
diameter from 10 to 50 microns, an increase of 0.615 in the 
value of Em. As stated previously, for c = 8 feet {xp=\790) 
in figure 7 (c), the increase in the value of Em is 0.170 for a 
corresponding increase in the value of the free-stream 
Mach number from 1.1 to 2.0. For a constant value of 
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(b) Effect of thickness ratio. Droplet diameter, 20 microns; free-stream Mach number, 1.4; pressure altitude, 15,000 feet; free-stream droplet 

Reynolds number, 453. 

(c) Effect of free-stream Mach number. Droplet diameter, 20 microns; airfoil thickness ratio, 0.02; pressure altitude, 15,000 feet. 

Figure 7. — Continued. Total collection efficiency of symmetrical double-wedge airfoils as function of scale parameter. Free-stream static 

temperature, 440° R; angle of attack, 0°. 
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(d) Effect of droplet size and chord length. Free-stream Mach number, 1.4; airfoil thickness ratio, 0.02; pressure altitude, 15,000 feet. 

(e) Effect of droplet size. Free-stream Mach number, 1.3; airfoil thickness ratio, 0.06; pressure altitude, 15,000 feet. 

Figi^re 7. — Concluded. Total collection efficiency of symmetrical double-wedge airfoils as function of scale parameter. Free-stream static 

temperature, 440° R; angle of attack, 0°. 
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cliord lengtli, varying the pressure altitude lias a relatively 
small effect on the total collection efficiency (fig. 7 (a)). 

Total collection efficiency as function of angle of attack. 
The previous discussion is concerned with the total col- 
lection efficiencies for the doulile- wedge airfoil at zero 
angle of attack only. The effect of angle of attack on the 
total collection efficiency is presented in figure 8 for aii-foil 
thickness ratios from 0,01 to 0.08, droplet diameter of 20 
microns, free-stream Mach number of 1.4, pressure altitude 
of 15,000 feet, and chord lengths of 1, 2, 4, and 12 feet. 
The range of angle of attack is from zero to tan“^(T/c). 

The total collection efficiency decreases slighth^ as the 
angle of attack increases. The slope of the curve of 
against tan a is zero at a = Q, by virtue of the symmetry of 
the double-wedge airfoil. As the angle of attack increases, 
the slope of the against tan a curve becomes negative, 
the rate of change becoming more pronounced as the angle 
of attack increases. The effect of increasing the chord length 
of the double-wedge airfoil is to decrease slightly the rate of 
decrease of Em with respect to tan a. A significant point 
that is well-illustrated in figure 8 is that, for a constant value 
of chord length, there apparentl}^ exists an envelope of the 
family of Em against tan a. curves that have the airfoil thick- 
ness ratio Tjc as the parameter. This envelope curve 
presents, in terms of the angle of attack, the maximum total 
collection efficiency that can be obtained for a symmetrical 
double-wedge airfoil of any thickness ratio, where the 
droplet size, free-stream Mach number, pressure altitude, 
and chord length are considered to be constant. 

For the double-wedge airfoil at a free-stream Mach 
number greater than the shock-wave-attachment Mach 
number, the decrease in total collection efficiency with an 
increase in angle of attack is o])posed to the trend experienced 
])y rounded-leading-edge airfoils at subsonic Mach numbers 
(irrespective of the S3nnmetry of the airfoil). For the latter 
t\^pe of airfoil at subsonic speeds, the increase in total 
collection efficiency (as defined herein) with increasing angle 
of attack is accounted for by the greatl\^ increased impinge- 
ment that occurs on the lower surfaces of these airfoils at 
angle of attack. The reduction of inij)ingement occurring 
on the upper surfaces of these airfoils at angle of attack is 
more than balanced by the increased ini])ingement occurring 
on the lower surfaces. On the other hand, for a given double- 
wedge airfoil at supersonic speeds with an attached shock 
wave, the rate of increase with angle of attack of the impinge- 
ment on the lower surface is less in magnitude than the rate 
of decrease with, angle of attack of the impingement on the 
up])er surface. This general trend for the double-wedge 
airfoil can be explained with the aid of figure 5 (a), in which 
the effect of tan a (tangent of the seniiapex angle) on the 
.total impingement rate is presented for a wedge of infinite 
extent. 

The forward upper surface of a s> mmetrical double-wedge 
airfoil effectively acts as the finite portion of one surface of an 
infinite wedge that is decreasing its semiapex angle (effect- 
ively" the thickness of wedge) as the angle of attack of the 
double-wedge airfoil increases. The forward! lower sui-face 
of this type airfoil effectively acts as the finite portion of one 
surface of an infinite wedge that is increasing its semiapex 



Figure 8. — Variation of total collection efficiency of symmetrical 
double-wedge airfoils with angle of attack. Free-stream static 
temperature, 440° R; droplet diameter, 20 microns; free-stream Mach 
number, 1.4; pressure altitude, 15,000 feet. 

angle as the angle of attack increases. In figure 5 (a) the 
increase in which is exactly ])roportional to the total 
impingement rate, becomes smaller with an increase in tan a 
(i.e., d^^mld{tan o-)^<^0) for all possible values of the semiapex 
angle a. F or example, in figure 5 (a) for the curve of Mi=l .4, 
the impingement on a sy^mmetrical double-wedge airfoil of 
0.08 thickness ratio and large chord extent at zero angle of 
attack would be very nearly proportional to twice the value 
of ^m at tan a of 0.08 (2f„j=0.146 ft). The impingement for 
the same airfoil under the same conditions and at an angle of 
attack of 2.291° (tan a=0.04) would be very nearly pro- 
portional to the sum of ^m at tan a of 0.04 and 0.12 
(fm.« + fm.i = 0.046 + 0.0905 = 0.1365 ft). 

Total collection efficiency as function of relative modified 
inertia parameter.— For air-flow fields that contain a shock 
wave, such as those considered herein, the water droplets 
upon crossing the shock wave suddenly have a velocity 
relative to the air. For the double-wedge airfoil with an 
attached shock wave, this initial relative velocity is the same 
for all droplets entering the air-flow field downstream of the 
shock wave. As discussed in previous sections, this common 
initial velocity may be considered the initial velocity of 
projection of a droplet in a reference frame having no air 
motion in it and moving relative to the fixed fiame of refer- 
ence at the constant velocity of the air downstream of the 
shock wave. Therefore, it is of interest to define a set of 
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inertia ])arameters, F and /'o, based on tlie motion of tlie 
dro])let in this moving reference frame. The relative 
inertia parameter based on the relative velocity Ui is defined 
as 

^ i ^ Pw^ F i ^ m . s 

~ fjLiC ~ c 

^vhere K is the inertia ])aranieter in the fixed frame of refer- 
ence, defined as 


Vi 

9 piC ~ 


(24) 


and x„i^s is that value of the maximum distance of travel 
obtained when Stokes’ law is assumed for the drag force on 
the droplet. 

The relative modified inertia ])arameter, also based on the 
relative velocity f/<, is defined as 


where Xf,i is defined in equation (9), and experimental values 
are used for the drag force on the dro])let. 

For the ])roblem of water-dro])let imjiingement on airfoils 
having rounded leading edges at subsonic s])eeds, it is ini- 
])ossible to define a similar relative modified inertia param- 
eter. Also, for any airfoil at siqiersonic speeds, it is impos- 
sible to define a relative modified inertia parameter, unless 
the shock wave from the leading edge is of constant strength 


in the vicinity of the airfoil, as in tlie case of the double- 
wedge airfoils with attached shock waves. If the strength 
of the shock wave is not constant, as for the case of a bow 
shock wave preceding a blunt obstacle, then the value of Fq 
is not constant for all the trajectories of a given dro])let size. 
Thus, the problem of impingement on wedges and double- 
wedge airfoils at supersonic speeds with attached shock 
waves is unique in that a relative modified inertia parameter 
can be defined and used in correlating the total collection 
efficiencies. 

The correlation of the total collection efficiency with the 
1 ‘elative modified inertia parameter Fq is shown in figure 9 (a) 
for di*oplet diameters varying from 2 to 100 microns and 
for pressure altitudes of sea level, 15,000, and 30,000 feet. 
For the symmetrical double-wedge airfoil of 0.06 thickness 
ratio at d/i = 1.3 and zero angle of attack, the droplet free- 
stream Keynolds number Rei varies from 42.1 to 2104. The 
values of for a given thickness ratio, Mach number, and 
angle of attack generally form the basis for a single curve 
with a small amount of scatter existing in the higher range 
of value of Fq. In the lower range of /'o, for all values of /?Ci 
the ])lotted ])oints have negligible scatter. The small scatter 
observed is iiossibly due to the existence of a very slight 
trend of the curves of E„, against Fq with Re^ for the higher 
values of Fq. 

The effect of increasing the thickness ratio is presented in 
figure 9 (b) for symmetrical double-wedge airfoils of 0.02, 
0.06, and 0.12 thickness ratio. Increasing the thickness 
ratio displaces the against Fq curve toward larger values 



(a) Correlation of pressure altitude and droplet size. Free-stream Mach number, 1.3; airfoil thickness ratio, O.Oh. 

Figure 9. — Total collection efficiency of sjunmetrical double-wedge airfoils as function of relative modified inertia parameter. Free-stream static 

temperature, 440° R; angle of attack, 0°. 
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(b) Effect of thickness ratio. Droplet diameter, 20 microns; free-stream Mach number, 1.4; pressure altitude, 15,000 feet. 

(c) Effect of free-stream Mach number. Droplet diameter, 20 microns; airfoil thickness ratio, 0.02; pre.ssure altitude, 15,000 feet. 

Figure 9. — ('oncluded. Total collection efficiency of symmetrical double-wedge airfoils as function of relative modified inertia parameter. Free- 

stream static temperature, 440° R; angle of attack, 0°, 
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of /^’o- (Miangin^ the thickness ratio of the double-wed^e 
airfoil does not alter the sha])e of the curve itself. 

The elfect of the free-streain Mach niiinher on the variation 
of K,n with /’(, is ])resented in figure 0 (c) for the symmetrical 
douhle-wedge airfoil of 0.02 thickness i*atio for dro])let diam- 
eter of 20 microns, free-stream Macli numbers of 1.1, 1.2, 
1.5, and 2.0, and ])ressure altitude of 15,000 feet. Increasing 
the Mach number dis])laces the entire curve of E,n against 
Fq toward smaller values of Fq. As the value of d/i increases, 
the I’ate of dis])lacement of the curve with increasing JA 
becomes smaller. The dis])lacement of the curve obtained 
1)\^ increasing the Mach numher from 1.1 to 1.2 is moi’e than 
that obtained b\' increasing it from 1.2 to 1.5 and from 1.5 
to 2.0. 

The relative inertia ])arameter /^Aind the relative modihed 
inertia ])arameter F^ of the moving frame of reference cor- 
res])ond, resiiectivel}’^, to the inertia jiarameter K and a 
modified inertia ])arameter Kq of the fixed frame of reference. 
The modified inertia ])arameter is defined in reference 14 as 

A'o=^ K (25) 

whei’e X is the maximum distance of travel of a dro])let 
])rojected into still air with the free-stream velocity IV 
The term X., is the value of the maximum distance of travel X 
when Stokes’ law is assumed for the drag force on the dro])let. 
The total collection efliciency of a symmetrical double-wedge 
airfoil is ])resented in a])])endix I) as a function of the modified 
inertia ])arameter /vq. 

SUMMARY OF RESULTS 

This re])ort ])resents an analysis of the })roblem of im- 
pingement of water drojilcts on a wedge and a do\d)le-wedge 
airfoil at suj)ersonic sj)eeds with attached shock waves. 
A\li(‘n a suitable em])ii‘ical relation is used for the drag co- 
edicient of a sphere, the analysis allows a closed-form inte- 
g]*ation of the ecpiations of motion for the water (Iroi)lets. 
The integration results in anal 3 dical expressions for the 
ecpiation of the trajectories, the dro])let velocit\^ at imy ])oint 
on the trajectories, the local impingement efliciencies, and 
the total rate of im])ingement. The results of the calcula- 
tions of rate, extent, and distribution of infringement of 
water dro])lets on wedges and symmetrical double-wedge 
airfoils are summarized briefh' as follows (tlie Mach number 
referred to is the free-stream Mach number, wiiich is greater 
than the attachment MacJi number for the wedge or the 
double-wedge airfoil) : 

1. At a given value of Mach number, dro])let size, and 
jiressure altitude, the local impingement efliciency as a 
function of the dimensional surface distance is the same for 
both the wedge and the s^mimetrical do\d)le-wedge airfoil 
at zero angle of attack, provided the tangent of the semiapex 
angle of the wedge is e([ual to the thickness i*atio of the sym- 
metrical double-wedge airfoil. 


2. For any Mach number, pressure altitude, and droplet 
diameter, the value of is equal to the sine of the semiapex 
angle for wedge or symmetrical double-wedge aii-foil. [^q 
is the value (maximum) of local impingement efliciency 
as distance from leading edge to ])oint of impingement 
a])i)roaches zero.] 

3. The effect of tiie free-stream static temperature on the 
local impingement efficiency and total collection efliciency 
is negligible for teni])eratures from 420° to 400° R. 

4. At constant values of Mach number, di-oplet size, and 
semia])ex angle of tlie wedge or corresponding tliickness 
ratio of the s^mimetrical double-wedge aii’foil, an increase in 
the pressure altitude increases slightly the local impingement 
efficiencies and total collection rates on wedges and sym- 
metrical double-wedge airfoils, but has a negligible effect 
on the surface extent of perceptible impingement. 

5. At constant values of Mach number, dro])let size, and 
pressure altitude, increasing the thickness ratio of the sym- 
metrical double-wedge airfoil or corresponding seniiapex 
angle of the wedge increases tlie local inii)ingement efficiency, 
has a negligible effect on tlie surface extent of jierceptible 
impingement, and decreases tlie total collection efficiency 
of the symmetrical double-wedge airfoil. 

0. At constant values of droplet size, pressure altitude, 
and semiapex angle of the wedge or thickness ratio of the 
s^unmetrical double-wedge airfoil, an increase in Mach 
number increases both the surface extent of impingement 
and the value of the local impingement efficiency, 

7. At constant values of jiressure altitude, seniiajjex 
angle of the wedge or thickness ratio of the symmetrical 
double-wedge airfoil, and Afach number, an increase in 
the (Iroj)let size increases considerably the surface extent of 
I)erceptible impingement, the local imjiingement efficiency, 
and the total impingement rate. 

8. The variation of total collection efficiency of the s\mi- 
metrical double-wedge aii-foil at zero angle of attack as a 
function of the scale parameter for constant values of tlie 
(Iroj)let fi’ee-stream lAwnolds number is similar in foi-m to 
that for subsonic airfoils. 

9. The total collection efficiency of the symmetrical 
double-wedge airfoil decrease's slightly as the angle of attack 
increases. 

10. For a symmetrical double-wedge airfoil of a given 
thickness ratio and Macli number, the values of total col- 
lection efficiency for a wide range of values of droplet free- 
stream Reynolds number comprise a single curve when 
plotted against the relative modified inertia parameter. Tlie 
effect of increasing tlie tliickness ratio or decreasing the 
Macli number is to displace the entire curve in the direction 
of larger values of relative modified inertia ])arameter. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautk’s 
Cleveland, Ohio, April 21, 
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SYMBOLS 

following symbols are used in this ro])ort: 7 

droplet radius, ft e 

drag coefficient, dimensionless 

airfoil chord length, ft f 

drag force on spherical water dro])let, lb 
droplet diameter, microns 

total collection efficiency (defined by eq. (22)), 
dimensionless 


^^3 


U 

U 

ir 

w 


/3o 


1 j* • i. • 4. 2 Ptod i ^m,s 

relative inertia parameter, - — 


9 PiC 


limen- 


sionless 


relative modified inertia ])arameter, — - 

m 

acceleration due to gravity, 32.2 ft/sec^ 


) ( 


^ m 
C 


inertia parameter, — 


2 pwd^Vi X 


PiC 


=—) dimensionless 
c 


modified inertia parameter, :^7^=-Kliniensionless 

X5 c 

Mach number 

empirical constant (used in eq. (5)), 2/3 
constants of flow field (defined by ecjs. (ISa), (bSb), 
and (18c), respectively) 
static pressure, Ib/sq ft 
gas constant, 53.3 ft-lb/(lb)(°F) 
droplet Reynolds numlier relative to air behind 
shock wave, 2 ap 2 Ulp 2 , dimensionless 
free-stream droplet Reynolds number, 2apiTq/Mi> 
dimensionless 

distance to point of impingement measured from 
leading edge for water dro])let that enters flow 
field behind shock wave at distance f above 
leading edge (eq. (15) and fig. 1), ft 
maximum thickness of airfoil, ft 
time, sec 

free-stream static temperature, °I{ 
magnitude of droplet velocity relative to air 
velocity downstream of shock wave, | Tsl, f])s 
velocity, fps or mph 
rate of water catch, lb/(hr)(ft span) 
local rate of water catch, ll)/(hr)(sq ft surface) 
liquid-water content of air, g/cu m 

J" U (It, displacement of droplet in moving frame 

of reference (relative to air stream); where 
the instant water droplet intercepts shock wave 
maximum value of displacement x as r^cx> as limit 
(maximum distance of travel of dro])let projected 
into still air with relative velocity f/d, ft 
value of maximum distance of travel obtained 
by assuming Stokes' law for drag force on 
droplet, ft 

angle of attack of airfoil, deg 
local impingement efficiency, d'^jdS, dimensionless 
maximum value of ^ as i0o = sin a, dimensions 
less 




P 

Pw 

(T 


4^ 

n 


ratio of specific heats (1.4 for air) 
empirical constant (used in eq. (5)), 0.158, dimen- 
sionless 

initial displacement of droplet trajectories from 
leading edge in direction normal to free-stream 
direction (eq. (13)) 

maximum value of initial droplet displacement f 
obtained when trajectory is tangent to wedge 
surface (theoretically as S^oo) 
absolute values of initial displacement from leading 
edge (in direction normal to free-stream direc- 
tion) of droplet trajectories that impinge at 
shoidder of upper and lower surfaces, respec- 
tively, of double-wedge airfoil 
distance along shock wave measured from wedge 
apex to point where droplet trajectory intercepts 
shock wave 
shock-wave angle 

angle between free-stream velocity vector and 
droplet velo(*ity vector Va 
maximum distance of travel of droplet projected 
into still air with free-stream velocity V,, ft 
value of distance of travel X obtained by assuming 
Stokes' law for drag force on water droplet, ft 
dynamic viscosit}^ based on static temperature, 
(lb)(sec)/sq ft 

angle between free-stream velocity vector V, and 
initial relative velocity vector (defined by eq. 
(13a)) 

components of droplet displacement referred to 
wedge surface (defined by eq. (12)), ft 
mass density of air, slugs/cu ft 
water-droplet mass density, 1.9398 slugs/cu ft 
semia])ex angle 

angle between dro])let velocity vector Va and 
velocity vector V 2 

dimensionless time variable, (S P 2 I Pud^)i 
phase angle, tan“^(/?Cr, 


r^/3e-^/^), 0^^g|(eq. (9)) 


scale parameter, (9cla)(pi/puj) 

ratio of air velocity downstream of shock wave to 
free-stream velocity V 2 IV, 

CO ratio of initial droplet relative velocity' to free- 

stream velocity UijV, 

vSubscripts: 
d dro])let 

i initial (at shock wave) 

im impingement 

I lower 

m maximum 

u upper 

1 free stream 

2 downstream of shock wave 

reference) 

Barred symbols denote vectorial quantities 


(in fixed frame of 
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APPENDIX B 


CLOSED-IORM INTEGRATION OF EIJUATION OF MOTION (RELATIVE TO AIR FLOW DOWNSTREAM OF SHOCK WAVE) 

'riu' vaviovis sU'ps lu'ccssury lo the closed-form integration 
of (‘<|ualion ((i), which is llie didVrontial oqualion of niotioii 
of the water droplet in ttie moving referenee frame, are given 


herein. 'Pl\e differential erpiation of motion ((>) can be re- 
written as 


(11 TI 2 24f , . Amr^rn vi\ 

. 11 — A. «■ <'+<•>'' « > 

(Bl) 

wluTo, for (■oiivoiiicnee. 


Rer=A,ra 

(B2) 


(B3) 

M2 


4 '3 p2 

-^12 — 77 

(B4) 

8 Pw 


Algebraic simiilincation and rewriting ri'sult in 

tlio form 

,/ -I 3 

T(i 

(B5) 

wlu're 


4 Ai 

A 2 24 

(Bb) 

and 


4 i m 7H 

Ai=eAiG 

(B7) 


E(|ualion (Bo) is not readily integrable in its j)resent form. 

J r ni 

]) = L 


JW lotting 


ami 


j T m 

I’H-r- 


o(|uation (Bo) boooinos 


^ <lp'^ A^Pi dp' 

' m p' m \+A^pip' 


(BS) 

(BO) 

(BIO) 


This form of the etpiation is readily integrated, and upon 
resubstiliiting the relations (BS) and (B9), there is obtained 
the following e.\j)ression for the velocity of tlie droplet as a 
function of the time (; 


r= r'ijexp j^(/— /fi) •"'-^‘'1 


1 / m 


(BID 


wlioro 7b is tho constant of integration. 
(lx 


Since r = (Miuation (Bit) (‘an be rewritten in integral 


form as 


JVt jexp[(/-/A) 


-\Jm 


dt + Bo (B 12 ) 


where 7b is the s(H‘oiuI integration constant. Ecination 
(Jil2) re([uir(‘s simjilifieation before a closed-form int(‘gration 
can be iierformed. Consider tlie following substitutions: 


m 

- 13 


(B13) 


.u=— ei(f 7 A 4 )-'^”'=— (Ad-"'" 

m ni ^ 


Ae=(rrA,)-'ex|) [-/bw/Aal 


(B14) 

(Bl5) 


By using the substitutions given by (Bl3), (Bl4), and (Bl5), 
it is i)ossible to wi‘ite equation (Bl2) as 


J r/x = Ja5 (.16C‘'- l)-""V///-b Ik 


(Bits) 


Before further steps can be taken in the closed-form integra- 
tion of the eejuation of motion of the droplet, the value of 
the empirical constant m must be determined. It is noted 
tliat in a})})roximating the curve of the drag coefHcient as a 
function of the local Keymdds nuinlHu* by a iTlation of the 
form given by 

(h) 


C/) — (1 -f e7^c7) 


it is possil)le to consider that the value of m is 2/3 and the 
value of c is 0.158. That the approximation of the experi- 
mental curve by the empirical relation is v(u-y good in the 
range of Kiwnolds numbers from about 0.5 to 500 can be 
seen by referring to figure 2, which ])resents a graph of the 
em])irical relation along with the drag-coefficient data of 
references 4 and 12. 

The use of m = 2/3 in (»(j nation (BlO) along with the sub- 
stitution 

e^=q (Bl7) 

and the use of formulas of integration given on pages 10 and 
17 of reference 15 allow equation (BlOj to be integrated as 
shown in the following stej)s: 


J'"-/ 


\,dq 


1)3/2 


-Ik 




x~/k 

2As ^ 


1 

^ 1 


^Aeq-l (B18) 


Jf the original time or independent variable t is reintroduced 
into equation (BI8) and the substitutions are made for /I5 
and Aq, the e(|uation has the form 




(A, 7/A2/3)i/2 


3.I3 






2/3 


(Bit)) 
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Substituting for and A 4 (except in the exponent of e) in 
equation (Bl9) results in 


stream of the shock wave as follows (in the frame of refer 
ence moving with the constant velocity 1 ^ 2 ).' 



V 1 

1 I 

to 

1 ^ 

1 


! /fe,. ^5, “ V”'’ ! 


(B20) 


Equation (B20) is the integrated equation with undeter- 
mined integration constants for the motion of the watei* 
droplets relative to the air velocity behind the shock wave. 
The integration constants are determined from the boundary 
conditions, which are 


U=Ui) 

x=^0 ) 


at ^ = 0 


The substitution of the boundary conditions, and thus tlie 
determination of the integration constants 7A and B 2 , 
results in the final form of the integrated ecpiation of motion 
for the water droplets relative to the air velocity down- 


6 p2 L ’ 






(7) 


The final form of the corresponding equation for the relative 
velocity of (he water droplet as a function of the dimension- 
less time variable is obtained from equation (Bll) as 


^ = Ii~V2 l(IiCr.r^'^e-^+ (8) 

It can be seen from equations (7) and (8) that as t approaches 
infinity the value of Cf approaches zero and that a limit 
exists for the value of x as ^ approaches infinity. This 
limiting value of x is 

^ (9) 

where 


APPENDIX C 

COMPARISON OF COLLECTION EFFICIENCY AND IMPINGEMENT RATE FOR SYMMETRICAL DOUBLE-WEDGE AIRFOIL 

WITH THOSE FOR NACA 0006^64 AIRFOIL 


A comparison of the total collection efficiency as a func- 
tion of a modified inertia parameter /C) for a symmetrical 
double-wedge airfoil at siqiersonic speeds (attached shock 
wave) with that for the NACA 000()-()4 airfoil (ref. 9) at 
free-stream Mach numbers less than critical is presented in 
figure 10. Both airfoils have O.OO thickness ratios and are 
at zero angle of attack. It is important to kee]) in mind 
that figure 10 does not allow a coni])arison of the two air- 
foils at the same Mach number, and such a comparison caii- 



Figure 10. — Comparison of total collection efficiency as function of 
modified inertia parameter for symmetrical double-wedge airfoil at 
supersonic speeds with attached shock wave with that for NACA 
0006-04 airfoil at free-stream Mach number less than critical. Air- 
foil thickness ratio, 0.06; angle of attack, 0°. 


not be made, because the analysis for the NACA 0006-64 
airfoil is not valid above the critical Mach number and that 
for the symmetrical double-wedge airfoil is not valid below 
the shock-wave-attachment Mach number. It is, however, 
of value to consider a conqiarison of the two airfoils with 
each operating within its ap])ro])riate speed range. In figure 
10 the rate of increase in with respect to is greater 
for the symmetrical double-wedge airfoil than for the NACA 
0006-64. For a constant droplet size, pressure altitude and 
temperature, and free-stream Mach numlier, resulting in a 
constant value of X, Kq varies inversely as the chord length 
{Kq=\/c). In general, therefore, the rate of decrease in 
with increasing chord length is greater in magnitude for the 
symmetrical double-wedge airfoil than for the NACA 0006-64 
airfoil. 

A comparison of the rate of total water catch per unit 
span for the symmetrical double-wedge airfoil with that for 
the NACA 0006-64 airfoil can be obtained if values (neces- 
sarily different) for the free-stream Mach numbers for the 
two airfoils are assigned. Assume that both airfoils are of 
the same thickness ratio and chord length and that they 
encounter identical icing conditions of droplet diameter (20 
microns), pressure altitude and temperature (15,000 ft and 
440° R), and liquid-water content (0.5 g/cu m), but have 
free-stream Mach numbers of 1.5 and 0.75, respectively. 
For the static temperature assumed, the Mach numbers of 
1.5 and 0.75 correspond to speeds of 1051 and 526 miles 
])er hour, respectively. Therefore, the magnitudes of the 
inertia parameter K and the free-stream droplet Reynolds 
number Rci for the double-wedge airfoil are twice as great 


Total collection efficiency, ^ Total collection efficiency, 


22 


REPORT 1159 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 





(a) Effect of j)ressure altitude and droplet size. Free-streaiii Mach numl)er, airfoil thickness ratio, 0.00. 

(1)) Correlation of thickness ratio. Droplet diameter, 20 microns; free-str(‘am Mach number. 1.4; i)ressure altitude, 15,000 feet. 

(c) Effect of free-stream Mach number. Droplet diameter, 20 microns; airfoil thickness ratio, 0.02; pressure altitude, 15,000 feet. 

Fir, THE 11. — Total collection efficiency of symmetrical double-wedge airfoils as function of modified inertia parameter. Frec-streani static 

temperature, 440° R; angle of attack, 0°. 
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as for the NACA 0006-64 airfoil. Varying the chord length 
fi'oin 1 to 20 feet produces a change in the value of Kq from 
0.386 to 0.0193 for the double-wedge airfoil and from 0.266 
to 0.0130 for the NACA 0006-64 airfoil. This variation in 
Ko for both airfoils results in values of of the same order 
of magnitude. P^or the given icing conditions, the following 
table lists for chord lengths of 4 and 20 feet the various per- 
tinent parameters and variables, including the rate of total 
water catch on the airfoil per unit span for both airfoils: 



Symmetrical double-wedge airfoil at A/i = l.5 (1 

’i = 1051 mph) 


Chord 
length, 
c, ft 

Free-stream 

droplet 

Reynolds 

number, 

Re, 

(a) 

Inertia 

parameter, 

K 

Distance 

ratio, 

X/X. 

(b) 

Modified 

inertia 

parameter, 

Ao 

(b) 

Total col- 
lection 
efficiency, 
Em 

Rate of 
water 
catch, ir, 
lb/(hr) 
(unit span) 

4 

20 

485 

485 

0. 542 
.108 

0. 178 
.178 

0. 0005 
.0193 

0.432 
. Ill 

17.9 

23.0 


NAC.\ 0000-04 airfoil at Mi =0.75 (Vi = 

520 mph) 


4 

243 

0. 271 

0. 240 

0. 0005 

0. 330 

0. 85 

2t) 

243 

.0542 

.240 

.0130 

. 128 

13.3 


“ Values of viscosity of air obtained from linear variation with temperature given in ref. 10. 
b Calculated as indicated in appendi.x D. 


The rate of total water catch on the airfoil per unit span is 
calculated from 

where T is in feet, w in grams per cubic meter, and Vi in 
miles per hour. P"or this particular example, the table shows 
that for the 4-foot chord, the double-wedge airfoil has a 
somewhat larger value of E,n than the NACA 0006-64 air- 
foil; and for the 20-foot chord, the opposite is the case. 
Comparison of the two airfoils for a given chord length 
shows that, as expected, the effect of the total collection 
efficiency on the rate of water catch is small when compared 
with the free-stream velocity ratio chosen (2:1). The most 
significant comparison to be obtained from the table is that, 
for the double-wedge airfoil, increasing the chord length by a 
factor of 5 (from 4 to 20 ft) results in an increase of only 28 
])ercent in the rate of water catch; whereas, for the NACA 
airfoil, a like change in the chord length results in an increase 
of 94 percent. This difference is the result of the fact that, 
unlike the local impingement on the rounded-leading-edge 
airfoils at free-stream Mach numbers less than critical, the 
local impingement at a given point on the surface of a 
double-wedge airfoil does not vary with the chord length of 
tlie airfoil, which is at a supersonic Macli number above the 
shock-wave-attachment Mach number. 


APPENDIX D 


TOTAL COLLECTION EFFICIENCY OF SYMMETRICAL DOUBLE- WEDGE AIRFOILS AS FUNCTION OF MODIFIED INERTIA 

PARAMETER 


The total collection efficiency E„, is presented in figure 11 
as a function of the modified inertia parameter /lo, as sug- 
gested in reference 14. The relation for Ko is 

Ko=^K (25) 


where K is the inertia parameter, defined as 


..2 

9 nic 


(24) 


Since K is equal to X^/c (ref. 4), Kq may also be written as 



The free-stream droplet Reynolds numbers used herein in 
calculating Ko range up to 2104. The empirical drag law 
(eq. (5)) used in this report for the droplet motion in the 
moving reference frame is valid up to the Reynolds number of 
approximately 500. Therefore, for the values of Ko pre- 
sented in figure 11 (and also fig. 10), the values of X/X, were 
obtained from table 1 of reference 4. 

The curves of Em against Ko in figure 11 have exactly the 
same form as the curves of Em against /^o, except that the use 
of Ko results in a displacement of the entire curve toward 
larger abscissa values. The displacement of the Em against 
Ko curves obtained by varying the thickness ratio (fig. 1 1 (b)) 
or the free-stream Mach number (fig. 11 (c)) is less than 
the displacement of the Em against Fq curves obtained by the 
same procedures (figs. 9 (b) and (c)). 
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